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The inclusive nonleptonic transition rate of B meson in 
gluonic penguin channel is calculated in QCD using the heavy 
quark mass expansion. We found the branching ratio to be 
approximately 0.25%. 
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I. INTRODUCTION 

Recently, there has been a significant interest in the 
physics of beauty hadrons. From the theoretical view- 
point, its attractiveness is due to the fact that both 
perturbative and nonperturbative corrections are under 
control on a scale of the heavy b-quark mass. On the 
other hand, it is the B-physics where the Standard Model 
can be tested with a high precision and a possible New 
Physics is expected to be eventually caught. If a discrep- 
ancy between the predictions of the Standard Model and 
the data is established, it would be necessary to invoke 
the New Physics for its resolution. 

The aim of the present paper is to calculate a con- 
tribution to the nonleptonic inclusive decay rate of the 
B-meson due to the giuon penguin channel (b— S'Sg).This 
study is motivated by the existence of the semileptonic 
BR problem 0]. Despite many theoretical efforts, this 
problem so far remains unsolved within the Standard 
Model with the 20 % gap between the data and theo- 
retical predictions. In attempts to its resolution, it is 
essential to take into account all relevant contributions 
which could enhance the BRgj. On the other hand, the 
gluon penguin contribution to the nonleptonic inclusive 
width has not been yet calculated. It is our purpose to fill 
this gap in the literature and calculate the corresponding 
partial width beyond the perturbation theory. We find 
that no considerable enhancement of BRsj occurs due to 
this decay channel, in agreement with the common belief. 
Thus, if one is expecting some significant enhancement 
of the BRgj in the gluon penguin then an appearance of 
a New Physics is anticipated Q . 

To calculate the inclusive transition rate in the b— >sg 
channel nonperturbatively we use a powerful method of 



pre-asymptotic expansion It is based upon an in- 

corporation of elements of the heavy quark effective the- 
ory (HQET) into the Wilson operator product expan- 
sion (OPE). In this method, one uses the optical theorem 
to connect the decay rate P with the transition operator 
T defined by 

f{Q -^f^Q) = ij d^xT{H+f{x)H,ffm, (1.1) 

where H^ff denotes an effective weak Hamiltonian rel- 
evant for the Q f transition . The nonlocal opera- 
tor T is expanded in a series of local operators with the 
heavy quark mass serving as a large parameter of the ex- 
pansion. The long distance physics is parametrised by 
a set of matrix elements of the local operators over the 
hadronic states. The transition rate is calculated in a 
systematic way through an expansion in inverse powers 
of the heavy flavour quark mass. Thus, inclusive decay 
rates are determined in terms of this universal set of ma- 
trix elements without further phenomenological input. 

The structure of the paper is as follows. In Sec. II 
we present in short the method of pre-asymptotic heavy 
quark mass expansion. The inclusive width in the b— s-sg 
channel is calculated in Sec. III. The final Section IV 
provides a summary of our results. 



II. PRE-ASYMPTOTIC EXPANSION FOR 
INCLUSIVE B-DECAYS 

Th e starting point of our analysis is the expression 
(1.1) for the transition operator T. In the case b ^ sg 
decay the effective weak Hamiltonian is 
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Here, Gp is the Fermi constant, Vm and Vst stand for 
the CKM matrix elements. The effective three point in- 
teraction couples the right-handed b quark with the left- 
handed s quark and the gluon. The Wilson coefficient 
Cm/ is a result of integrati ng o ut hard gluons which gen- 
erate th e ef fective vertex (2.1). For the effective Hamil- 
tonian d^l) Cw{mt) ~ -0.150 The calculation of 
the decay rate P is based on the optical theorem, which 
relates the inclusive width to the imaginary part of the 
forward scattering amplitude: 
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The heavy meson mass Mb in Eq.(2^) accounts for the 
proper normaUsatfon. As was pointed out in Ref. 
the expression (2^) suggests a way for calculation of F 
through a systematic expansion in inverse powers of the 
heavy quark mass. The nonlocal operator T is expanded 
in a series of local operators, the operators of lowest di- 
mensions being Oi — bb and O2 — biaGb. Thus, T can 
be written in the form 



T ^ C\{B\Oi\B) + C2{B\02\B) 



(2.3) 



where dots denote contributions of higher order opera- 
tors. Since these terms are suppressed by powers of the 
heavy quark mass their contributions will be omitted in 
the following calculations. 

In the limit of heavy quark mass, the expectation value 
of operator Oi over a meson state is 
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The first term in the expression (2.4) is equal, up to nor- 
malisation, to the density of b quark in B meson. Tak- 
ing into account the normalisation of the B meson state 
{B\B) = 2Mb, we get 



{B\h^b\B) 
2Mb 
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The second and third terms in (2.4) are suppressed by 
inverse b quark mass squared. Note, that the second 
operator is just the O2 operator which appears in the ex- 
pansion for the Oi matrix element and plays the role of 
chromomagnetic i nter action The operator O2 arises 
in the expression (2^) for F also in a direct way through 
the expansion of the transition operator T. This hap- 
pens when the interaction of s quark and gluon with an 
external gluonic field is taken into account. The relevant 
propagators are 



S{x) = 
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Th e third operator on the right-hand side of equation 
( ^.4| ) describes the kinetic energy of the b quark inside 
the hadron. Its matrix element can be estimated to 
be 



{B\bn^b\B) 
2Mb 



0.5 GeV^ 



(2.9) 



III. GLUON PENGUIN CHANNEL IN 
INCLUSIVE DECAY 

In order to calculate the inclusive transition rate of B 



meson in the gluon penguin channel, we use Eq. (1.1) with 
the effective weak Hamiltonian defined in (p?l|). 

T = Ai I dxT{b{x) — ;r-^<Jal3Gaf:i(x)s{x) 



xs(0)a,,i^G,, (0)6(0)} 
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In the limit of the heavy quark mass, the b quark wave 
function is rescaled and the mechanical part is extracted 
Pi: 



b{x) = 6(0)e 



(3.2) 



Here p denotes the four-momentum of the b quark inside 
the hadron. The lowest order contributions to the ex- 



pression (3.1) for the transition operator T are given by 
the three diagrams (Figs. 0, ||). 
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FIG. 1. Zero order contribution to the transition operator 



T. 



Here G denotes the dual of the field strength. The ma- 
trix element of the chromomagnetic operator is propor- 
tional to the splitting between squared masses of B* and 
B mesons: 

^^|^-^(M|.-AfI).0.7Gey^ (2.8) 
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FIG. 2. First order contributions to the transition operator 
T: (a) the contribution of the quark propagator in an external 
field and (b) the contribution of the gluon propagator in an 
external field. 



IV. CONCLUSIONS 

In the present paper, we calculated nonperturbatively 
the inclusive nonleptonic transition rate of the B me- 
son in the gluonic penguin channel b — >sg by means of 
the pre-asymptotic heavy quark mass expansion. The 
obtained BR is about 0.25% while the nonperturbative 
contribution to the transition rate is about 4%. Our re- 
sult provides the Standard Model predictions beyond the 
perturbation theory and is to be checked experimentally. 
The number obtained means that a significant contribu- 
tion to the b — >sg decay mode is only possible within the 
New Physics scenarios. 
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These diagrams are obtained when the s quark and 
gluon fields are contracted in Eq. (3.1). Expressions (2.6), 



( p.7[ ) are used for the propagators in the external gluonic 
field. 

The first diagram (Fig. gives rise to the operator 
Oi while the other two (Fig. g|) produce the operator O2 ■ 
Th e ex plicit calculations for the coefficients Ci and C2 
in (2^) yield the following answers: 



C, = Ai — )CFmlH~ml) 
on 

^^=^[(2Ar)8^-8^]-^"(-'"^) 

Here, Nc stands for the number of colours and 

Cp = {N^-l)/2Nc. Using the fact that /mln(-m^) = tt 

we arrive at the final expression for the inclusive decay 

rate 
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Substituting expressions (^), (^, (|2^) and ( |2.9| ) 
for the matrix elements and using all numerical values 
we finally determine the transition rate: 



I.I X 10"^^ GeV. 



The appropriate branching ratio is 2.5x10 
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